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LONG-TERM  GOALS 

The  long-term  goal  of  “Imaging  System  Performance  and  Visibility  as  Affected  by  the  Physical 
Environment”  is  to  investigate  the  relationships  between  physical  forcing  and  optical  properties, 
including  the  imaging  performance  parameters.  Our  efforts  are  aimed  at  supporting  the  Radiance  in  a 
Dynamic  Ocean  (RaDyO)  program,  which  is  devoted  to  the  topic  of  light  propagation  and  imaging 
across  the  air-sea  interface  and  within  the  surface  boundary  layer  of  natural  water  bodies.  The  primary 
goals  of  the  RaDyO  program  are  to: 

1)  Examine  time-dependent  oceanic  radiance  distribution  in  relation  to  dynamic  surface  boundary 
layer  (SBL)  processes. 

2)  Construct  a  radiance-based  SBL  model. 

3)  Validate  the  model  with  field  observations. 

4)  Investigate  the  feasibility  of  inverting  the  model  to  yield  SBL  conditions. 

OBJECTIVES 

The  objectives  of  this  project  are  to: 

1)  Detennine  the  relationships  between  physical  processes  and  the  optical  properties  that  affect 
underwater  visibility  and  the  imaging  performance  parameters. 

2)  Year  2  Option:  Investigate  the  effects  of  vertical  variability  (e.g.,  vertical  layering)  of  optical 
properties  on  the  image  performance  parameters,  i.e.  modulation  transfer  function  (MTF)  and 
point  spread  function  (PSF). 


We  discuss  the  MTF  of  two  different  optical  water  types:  (1)  a  dynamic,  shallow-water,  eutrophic 
environment  and  (2)  a  relatively  clear,  deep  water,  mesotrophic  environment,  and  the  relationships 
between  the  MTF  and  optical  properties  and  physical  processes  (Chang  and  Twardowski,  2011). 
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APPROACH 


The  MTF,  which  is  the  magnitude  of  the  Fourier  transform  of  the  PSF,  was  computed  for  two  different 
water  bodies  using  measurements  of  optical  properties  and  analytical  formulations  described  by  Wells 
(1973)  and  Hou  et  al.  (2007).  The  MTF  of  an  aquatic  medium,  H(\|/,  z),  can  be  represented  by: 

H(\|/,  z)  =  e'D(¥)Z,  (1) 

where  Z  is  the  range,  D(\|/)  is  the  decay  transfer  function  (DTF),  and  \( /  is  spatial  frequency.  By 
application  of  small-angle  scattering  approximations  (Wells,  1973)  the  collimated  DTF  can  be 
approximated  as: 

Dc(¥)  =  c-S(¥),  (2) 

where  c  is  the  beam  attenuation  coefficient  and  S(i{/)  is  the  light  scattered  back  into  the  acceptance 
cone: 


0max 

S(¥)  =  2n  J0  Jo  (2tu0\|/)  P(0)  0  d0,  (3) 

where  (3(0)  is  the  volume  scattering  function  (VSF),  0  is  the  scattering  angle,  and  J0  is  the  zero-order 
Bessel  function.  The  desired  focused  DTF  can  then  be  obtained  by  transformation  in  terms  of  spatial 


frequency  (2tt0\|/)  (see  Wells,  1973  for  more  details): 

D(¥)=  C  Dc(\|/  t)  dt.  (4) 

In  situations  where  the  VSF  may  not  be  quantified,  the  following  assumptions  about  the  VSF  have 
been  made  (e.g.,  Wells,  1973): 

P(0)  =  (b  0O)  /  [2n  (  0O2  +  02)  3/2],  (5a) 

where  0O  is  mean  square  angle  and  b  is  the  scattering  coefficient: 

b  =  2n  jQ  P(0)  sin  0  d0;  (5b) 

therefore  a  closed-form  solution  for  the  MTF  is  possible  without  direct  measurements  of  the  VSF: 
D(¥)  =  c  -  [b  ( 1  -  e271 00  v)]  /  (2n  0O  ¥).  (6) 


The  optical  data  necessary  for  MTF  computations  were  collected  from  two  field  sites  as  part  of  the 
ONR  sponsored  RaDyO  program:  (1)  Scripps  Institution  of  Oceanography  (SIO)  Pier,  a  shallow-water, 
eutrophic  environment  and  (2)  the  Santa  Barbara  Channel  (SBC),  a  deeper,  mesotrophic  environment 
(Chang  et  al.,  2010).  The  data  processing  and  analysis  approach  utilized  factory-recommended  and 
published  calibration  and  correction  procedures  for  IOP  sensors  (e.g.,  Zaneveld  et  al.,  1994;  Pegau  et 
al.,  1997;  Sullivan  et  al.,  2006;  Zhang  et  al.,  2009).  Profiler  measurements  were  used  to  compute 
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optical  products:  backscattering  ratio,  bhp  (X)  =  bhp(A)/bp(A);  the  slope  of  the  cp(k)  spectrum,  y,  an 
indicator  of  particle  size  distribution  (Boss  et  ah,  2001);  the  real  part  of  the  index  of  refraction  of 
particles,  np,  from  y  and  bhp  (X)  (Twardowski  et  ah,  2001),  useful  for  obtaining  infonnation  about  the 

relative  density  of  particles;  the  single-scattering  albedo,  co0  =  b  /  c;  and  chlorophyll  concentration 
(Chi)  from  chlorophyll  fluorescence  and  spectral  absorption  using  the  methods  of  Sullivan  et  al. 

(2005).  The  MASCOT  data  were  used  to  compute  a  “bubble  index”  (BI),  where  BI  =  p(70°)  /  P(120°). 
This  BI  is  based  on  theoretically  computed  VSFs  of  relatively  large  bubbles  (>  10  pm),  which  show 
pronounced  enhancement  of  the  VSF  between  60°  and  80°  (Czerski  et  ah,  in  press;  Twardowski  et  ah, 
submitted  manuscript).  Hence,  increases  in  the  BI  can  be  indicative  of  a  relative  increase  in  the 
concentration  of  larger  bubbles. 

Relationships  between  physical,  hydrographic,  and  optical  properties  and  the  imaging  performance 
parameter,  MTF,  were  determined  for  the  SBC  experiment  by  use  of  wavelet  analysis.  Continuous 
wavelet  transforms  (CWTs),  cross  wavelet  transforms  (XWTs),  and  wavelet  coherence  (WTC)  are 
useful  tools  for  analyzing  localized,  intermittent  oscillations  in  data  series,  as  well  as  the  relationship 
and  linkages  between  two  data  series  and  the  strength  of  the  relationship(s)  (see  Grinsted  et  al.,  2004). 
Fourier  expansion  (e.g.,  fast  Fourier  transform;  FFT)  has  traditionally  been  the  technique  employed  for 
determining  the  frequencies  present  in  a  signal.  However,  FFTs  have  only  frequency  resolution  and  no 
time  resolution.  In  other  words,  although  FFTs  can  reveal  all  frequencies  present  in  a  signal,  they 
cannot  detennine  when  these  frequencies  are  present  in  time.  Similar  to  this,  coherence  analysis  can 
be  used  to  find  common  periodicities  between  two  different  time  series;  however,  determined 
commonality  is  not  localized  in  time  or  space.  Wavelet  analysis  can  be  utilized  to  expand  data  series 
into  time-frequency  space  and  thus  identify  localized  intermittent  periodicities  at  multiple  frequency 
and  temporal  resolutions. 

WORK  COMPLETED 

Measured  optical  properties  at  532  nm  from  SIO  Pier  and  the  SBC  were  used  to  compute  the  imaging 
performance  parameter,  MTF,  at  ranges  set  equal  to  the  water  depth.  The  MTF  fonnulation  presented 
by  Wells  (1973)  (Eqs.  1-4)  was  utilized,  with  near-forward  angle  scattering  provided  by  the  LISST- 
100X.  The  ^-approximation  (Eq.  5)  was  evaluated  using  SIO  Pier  data  by  comparing  the  MTF 
computed  using  the  LISST-measured  near-forward  angle  scattering  (9  between  0.05  and  8.25°)  and  the 
order  zero  Bessel  function,  J0  (Eq.  3)  versus  that  computed  using  an  approximated  VSF  (Eq.  5)  and  a 
mean  square  angle  determined  from  LISST  measurements  (0.07  rad). 

The  p -approximation  allowed  us  to  perform  a  sensitivity  analysis  of  particle  concentration  versus 
composition  effects  on  the  variability  of  the  MTF  in  order  to  determine  the  relative  importance  of 
particle  effects  on  the  MTF.  In  the  context  of  Eqs.  5  and  6,  the  particle  composition  effect  is  embodied 
in  the  mean  square  angle,  90.  Particle  concentration  effects  were  explored  by  varying  the  beam 
attenuation  coefficient,  cpg(532),  between  0.5  and  1 .0  nf1  by  steps  of  0. 1  m'1  while  keeping  the  mean 
square  angle,  0o,  constant  at  0.07  rad,  which  is  90  computed  for  the  LISST  measurements.  The  mean 
square  angle  was  then  varied  between  0.03  and  0.19  rad  by  steps  of  0.04  rad  while  keeping  cpg(532) 
constant  at  0.6  nf 1  to  investigate  relative  particle  composition  effects.  Note  that  we  are  using 
variability  in  90  as  a  proxy  for  relative  changes  in  particle  composition.  The  range  was  set  at  2  m  and 
co0  was  0.9  in  both  situations. 
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Relationships  between  the  MTF  and  other  parameters  were  quantified  using  a  series  of  Morlet  wavelets 
that  were  applied  as  band-pass  filters  to  2-m  time  series  of  the  MTF,  physical  properties,  and  optical 
properties  and  products.  The  wavelets  were  stretched  in  time  by  varying  their  scale  and  nonnalizing 
them  to  have  unit  energy.  The  statistical  significance  of  wavelet  power  was  assessed  for  various 
properties  to  detennine  dominant  scales  of  variability  and  the  periods  at  which  they  occur  (Grinsted  et 
ah,  2004).  Periods  and  frequencies  of  high  and  significant  wavelet  power  were  interpreted  as  dominant 
modes  of  variability  at  specific  times.  XWTs  and  WTCs  were  applied  to  time  series  data  to  expose 
high  common  power  and  relative  phase  between  two  variables  in  time-frequency  space.  A  Monte- 
Carlo  test  was  perfonned  to  detennine  the  statistical  significance  of  the  computed  coherence.  In 
addition  to  CWT,  XWT,  and  WTC  analysis  for  2-m  time  series  data,  WTC  methods  were  applied  to 
depth-resolved  profiles  of  physical  and  optical  properties  and  the  MTF.  Using  this  method,  common 
periodicities  between  two  different  parameters  in  vertical  space  were  evaluated  as  a  function  of  depth. 

RESULTS 

Variability  of  near-surface  MTF  at  both  field  sites  was  strongly  influenced  by  physical  forcing  - 
meteorological  processes  in  the  SBC  and  rip  currents  at  SIO  Pier.  The  2-week  SBC  field  experiment 
can  be  partitioned  into  three  distinct  periods:  (1)  diurnal  winds  (11-15  September),  (2)  persistent 
strong  winds  and  upper  water  column  mixing  (15-18  September),  and  (3)  an  advection  event 
characterized  by  relatively  high  Chi  and  elevated  optical  quantities  (19  -  21  September)  (Figure  1). 
During  diurnal  winds  and  the  sustained  strong  wind  event  in  the  SBC,  higher  wind  speeds  were 
associated  with  particle  characteristics  indicative  of  larger,  non-pigmented  particles  with  higher 
relative  indices  of  refraction  and  decreases  in  the  MTF.  The  presence  of  larger  particles  with  higher 
index  of  refraction  observed  during  the  high  winds  of  the  SBC  RaDyO  experiment  could  have  been  the 
result  of  wind-induced  mixing  of  the  upper  water  column  waters  and  surface  bubble  injection  or  an 
increase  in  non-chlorophyll  containing  particles  (perhaps  minerogenic  particles  or  senescent  diatoms) 
or  a  combination  of  both  of  these  factors.  Anderson  et  al.  (2008)  and  references  therein  have  reported 
on  phytoplankton  succession  in  the  SBC,  with  diatom  dominance  during  upwelling  periods  (spring  and 
summer),  followed  by  continued  surface  depletion  of  nutrients  and  stratification,  and  likely  diatom 
sinking  in  the  autumn.  It  is  possible  that  these  senescent  diatoms  were  retained  in  near  surface  waters 
during  strong  mixing  associated  with  high  winds. 
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Figure  1.  Time  series  of  (a)  wind  stress  (blue)  and  wind  speed  (red),  (b)  cpg(532)  (brown  dots)  and 
Chi  (green  triangles)  at  2  m,  (c)  near-surface  np  (blue  dots)  and  /(red  triangles),  near-surface  (d) 
volumetric  bubble  concentration  derived  from  BI,  and  (e)  MTF  for  W=  3  rad1  (black  dots)  and  *F= 
50  rad1  (gray  triangles)  calculated  from  2  m  optical  data  collected  in  the  SBC. 


Surface  bubble  injection  is  often  associated  with  strong  sustained  winds,  which  in  itself  would  result  in 
substantive  changes  in  particle  composition  characteristics  in  surface  waters  (Zhang  et  al.,  2011). 
Information  about  bubble  concentration  was  obtained  from  the  BI,  which  was  derived  from  the  ratio  of 
MASCOT-measured  VSFs  at  70°  and  120°.  The  BI  can  be  used  as  a  proxy  for  the  concentration  of  the 
large  bubble  subpopulation,  following  Twardowski  et  al.  (submitted  manuscript).  The  derived 
concentration  of  relatively  large  bubbles  (>  10  pm  in  diameter)  increased  by  nearly  an  order  of 
magnitude  on  15  September  2008  and  by  a  factor  of  three  on  18  September,  coinciding  with  periods  of 
sustained  elevated  wind  speeds  (Figure  Id).  The  increase  in  the  concentration  of  the  large  bubble 
subpopulation  was  accompanied  by  a  decrease  in  the  slope  of  the  particulate  attenuation  spectra, 
indicative  of  an  increase  in  particle  size  (Figure  lc). 

Although  not  statistically  significant  in  2-m  time  series  data,  the  observed  advection/bloom  event 
starting  on  19  September  2008  resulted  in  significant  coherence  between  the  MTF  and  the  IOPs  and 
IOP  products  indicative  of  particle  concentration  and  composition  at  about  20-m  water  depth.  This 
event  was  marked  by  a  decrease  in  the  MTF,  increases  in  the  magnitude  of  all  IOPs,  Chi,  and  np  and  a 
decrease  in  y,  indicative  of  an  increase  in  the  concentration  of  larger  biogenic  particles  with  relatively 
higher  refractive  index  (Figure  1).  The  concentration  of  the  large  bubble  subpopulation  during  this 
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time  period  was  relatively  low,  suggesting  lesser  effects  due  to  surface  bubble  injection  and  more  to  a 
phytoplankton  bloom  that  was  associated  with  a  shift  in  the  SBC  circulation  (Figure  1). 

Highly  reflective  minerogenic  particles  with  relatively  high  bulk  particle  index  of  refraction  were 
associated  with  SIO  Pier  rip  currents,  which  resulted  in  a  large  increase  in  the  MTF.  These  particles, 
although  not  explicitly  sampled  or  analyzed,  appeared  like  mica  flakes,  which  although  rare  in  beach 
sand,  have  been  observed  along  beaches  in  San  Diego,  CA.  This  increase  in  the  MTF  was  likely  due  to 
enhanced  reflectivity  coupled  with  relatively  high  visibility;  beam  attenuation  was  not  greatly  elevated 
during  the  rip  current  event  and  Chi  was  at  a  minimum  (Figure  2). 


Date  in  January  2008 


Figure  2.  SIO  Pier  time  series  of  2-m  (a)  cpg(532)  (brown  dots)  and  Chi  (green  triangles)  and  (b) 
bbp(532)  (blue  dots)  and  np  (red  triangles);  and  (c)  MTF  (W=  3  rad'1  in  black  dots 

and  tF=  50  rad1  in  gray  squares). 

These  observations  are  particularly  interesting  as  one  would  expect  the  imaging  perfonnance 
parameter  to  be  influenced  primarily  by  particle  concentration  as  revealed  by  the  magnitude  of  IOPs; 
however  we  found  that  particle  size  and  type  (relative  density  information  gained  from  the  real  index 
of  refraction  of  particles)  were  important  sources  of  variability  in  the  MTF. 

Our  results  show  that  the  effects  of  particle  composition  on  the  shape  of  the  VSF  in  the  near-forward 
direction  are  important  to  imaging  performance,  which  is  a  unique  observation  facilitated  by 
advancements  in  in  situ  optical  instrumentation.  By  performing  a  sensitivity  analysis  using  the  p- 
approximation  fonnulation  for  the  MTF,  we  were  able  to  elucidate  the  relative  importance  of  particle 
effects  on  the  MTF  (Figure  3).  Results  suggest  that  at  0O  between  0.07  and  0. 1 1  rad,  particle 
concentration  and  composition  affected  MTF  variability  similarly.  At  90  greater  than  0.1 1  rad,  MTF 
variability  was  stronger  with  variable  IOPs,  i.e.  particle  concentration  effects  dominated.  Strikingly, 
when  0O  is  between  0.03  and  0.07  rad  (i.e.  particle  size  distributions  trending  toward  larger  sizes),  the 
MTF  is  much  higher  for  frequencies  less  than  30  rad'1  (Figure  3).  Here,  particle  composition  effects 
dominated  MTF  variability,  in  a  relative  sense.  For  the  specific  ranges  of  variability  in  optical 
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properties  we  observed  during  the  study,  these  modeling  results  are  consistent  with  the  strong 
coherence  observed  between  particle  composition  parameters  and  imaging  performance. 
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Figure  3.  Particle  concentration  versus  composition  effects  on  MTF  variability,  (a)  MTF  computed 
by  varying  the  beam  attenuation  coefficient,  c  (as  indicated),  while  keeping  the  mean  square  angle, 
9„,  constant  at  0.07  rad.  (b)  MTF  computed  by  varying  6„  (as  indicated)  while  keeping  c  constant  at 
0.6  m1.  The  range  was  set  at  2  m  and  co,,  was  0.9  in  both  situations. 


The  small-angle  approximation  (Eq.  5;  e.g.,  Wells,  1973)  assumes  strong  forward  scattering,  which 
was  not  the  case  during  the  SIO  Pier  rip  currents.  To  explore  this  further,  we  evaluated  MTFs 
computed  from  the  DTF  using  measured  P  (Eq.  3)  versus  those  calculated  using  approximated  p  (Eq. 

6;  0O  =  0.07  rad)  for  the  SIO  Pier  data  set  (Figure  4).  The  p-approximation  method  produced  very 
similar  MTFs  to  those  computed  using  measured  P,  particularly  at  high  spatial  frequencies  (>  10  rad'1). 
MTFs  at  spatial  frequencies  less  than  10  rad'1  were  not  comparable  during  periods  when  the 
backscattering  coefficient  exceeded  0.012  m'1.  Results  from  the  two  methods  deviated  greatly  from 
each  other  at  all  spatial  frequencies  during  the  rip  current  event  (Figure  4).  Here,  we  show  that  in 
highly  backscattering  environments  such  as  those  with  substantial  minerogenic  particles,  the  shape  of 
the  VSF  cannot  be  assumed  when  computing  imaging  perfonnance  parameters  such  as  the  MTF. 

IMPACT/APPLICATIONS 

Optically-derived  particle  characteristics  such  as  bulk  particle  index  of  refraction  and  particle  size 
distribution  were  shown  to  be  significantly  related  to  the  variability  of  the  MTF.  This,  in  addition  to 
the  great  temporal  variability  in  environmental  conditions  observed  at  two  different  field  sites  stresses 
the  importance  of  obtaining  or  predicting  optical  properties  and  derived  particle  composition 
characteristics  (describing  relative  particle  size  and  type)  at  temporal  and  spatial  scales  necessary  to 
quantify  rapid  changes  in  image  processing  parameters.  Importantly,  estimates  of  water  clarity  or 
turbidity  are  not  the  only  factors  influencing  the  variability  of  the  image  performance  parameters, 
particularly  at  steep  near-forward  scattering  angles,  where  variability  in  the  MTF  was  found  to  be 
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dominated  by  particle  composition  effects.  We  also  show  that  the  standard  approximation  for  the  VSF 
cannot  be  used  to  compute  imaging  perfonnance  in  environments  with  highly  backscattering/reflective 
particles. 
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Figure  4.  SIO  Pier  time  series  of  MTF  at  *F=  (a)  3  rad1  and  (b)  50  rad1,  with  MTFs  computed 
using  the  /3-approximation  (e.g.,  Wells,  1973 )  shown  with  gray  triangles. 
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